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In view of recent experiments on ultra-cold polarized fermions, the zero-range potential approach is 
generalized to situations where two-body scattering is resonant in the p-wave channel. We introduce 
a modified scalar product which reveals a deep relation between the geometry of the Hilbert space 
and the interaction. This formulation is used to obtain a simple interpretation for the transfer rates 
between atomic and molecular states within a two branches picture of the many-body system close 
to resonance. At resonance, the energy of the dilute gas is found to vary linearly with density. 
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Studies of regimes with strong correlations is actu- 
ally one of the most challenging directions in the field 
of ultra-cold atoms. A spectacular example of such a sit- 
uation is given by the BCS-BEC crossover of a two spin- 
component Fermi gas, obtained via tuning the scattering 
length of binary collisions by using Feshbach resonances 
0. Another interesting example is resonant scattering 
in p-wave channels 0. [1 El 0. EJ with spin-polarized 6 Li 
and 40 K atoms which opens the possibility of a BCS-BEC 
crossover in higher angular momentum channels in the 
near future. The major interest in these systems stems 
from the fact that the strength of the correlations can be 
tuned arbitrarily while the mean inter-particle distance 
remains large with respect to the range of interatomic 
forces. Consequently, low energy properties are indepen- 
dent of the non-universal short range physics. Ultra-cold 
gases can therefore be considered as model systems for 
accurate studies of quantum many-body problems un- 
derlying many interesting phenomena in condensed mat- 
ter physics. For two spin-component fermions where the 
scattering amplitude in the s-wave channel dominates, 
several effective two-body potentials are used both in 
computational and analytical studies of the so-called uni- 
tary regime 0, IE IU El El 0], where the scattering 
length diverges. In this context, the zero-range pseudo- 
potential approach is very appealing as it captures all 
the effects of the interaction with only one parameter 
[Til fl3L IbH ] for broad s-wave Feshbach resonances, and 
with two parameters in narrow s-wave Feshbach reso- 
nances |l5| . Furthermore, this approach is particularly 
useful to obtain exact solutions of few body problems 
HHH [Hill, and also to improve unavoidable approx- 
imations in many-body systems [14L Il9| . 

In this letter, a general zero-range treatment of reso- 
nant p-wave interactions is introduced. This formalism is 
used in the context of ultra-cold spin-polarized fermions, 
where binary s-wave scattering is suppressed due to the 
Pauli principle. We show that the Hilbert space associ- 
ated with the zero-range Hamiltonian has to be defined 
with a new metrics. The modified scalar product illus- 
trates nicely an intrinsic connection between the geom- 
etry of the Hilbert space and the interaction. The for- 
malism is applied to an effective model to obtain some 



physical insight into the crossover regime of a many-body 
system made of particles interacting in pairs via the p- 
wave channel. A two branches picture is obtained for the 
equation of state in the neighborhood of the Feshbach 
resonance, giving a simple interpretation for the transfer 
rates between atomic and molecular states. At resonance, 
in the dilute regime, the model predicts a ground state 
energy which varies linearly with respect to the atomic 
density. 

Without any loss of generality, we introduce the for- 
malism by considering two identical particles in absence 
of an external potential. The particles of mass m, posi- 
tions (fi,r2) and relative coordinates f — ?\—?2 are de- 
scribed in their center of mass frame by the wave func- 
tion \?(r). The interaction between the particles occurs 
only in the p-wave channel of the two-body system and 
is modeled by the following zero-range pseudo-potential: 
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where 1Z[.] is a regularizing operator defined by: 
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In Eq.QJ, V s is the scattering volume. In the resonant 
regime, |V S | is arbitrarily large and in contrast with the 
analog unitary regime in the s-wave channel, a second 
parameter (denoted here by a) is essential for a descrip- 
tion of the shape of the two-body scattering amplitude 
pof and has been introduced in the expression of the 
regularizing operator. In Eq.Q, a surface integration 
is performed over the sphere S r of radius r centered at 
r = 0, d 2 fl is the elementary solid angle and e r = r/r. 
This integration ensures that the pseudo-potential acts 
only on the p-wave component of the wave function [2l| . 
As a consequence, the wave function has a 1 /r 2 singular- 
ity for r — y which is a general feature of the zero-range 
potential approach in the p-wave channel. The role of 
the pseudo-potential in the Schrodingcr equation is two- 
fold: first, it imposes a specific boundary condition for 
the wave function in the vicinity of the 1 /r 2 p-wave sin- 



2 



gularity: 
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and, second, it cancels a diverging term proportional to 
(V<5) coming from the action of the Laplacian on the 
wave function. Equation © is analogous to the contact 
condition introduced by H. Bethe and R. Peierls for a 
description of s-wave scattering |22f and represents an 
alternative way to formulate the zero-range approach. 
The pseudo-potential in Eq.ffl) generalizes other zero- 
range p-wave pseudo-potentials [23, |24( , where the reso- 
nant regime must be described through an energy depen- 
dent scattering volume. Contrary to these approaches, 
the pseudo-potential Q can be used directly in time de- 
pendent problems or in situations where the two-body 
collisional energy is not explicitly defined. 

As a first illustration of the formalism, we deduce the 
two-body eigenstates of energy E in their center of mass 
frame in absence of any external potential. For this pur- 
pose, we use the integral equation: 



*(f) = * (r) - Jd 3 r 'G E {r,r')(f' \V\V) , 
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where Ge is the one-body outgoing Green's function at 
energy E. For positive energies, E — h 2 k 2 /m, Eq.QJ 
yields the following expression for the scattering states: 



$ g (r ) = exp(ife.f) + 3V S e r .Tl[^ R ] d r 



exp(ifer) 
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Applying the regularizing operator TZ[ . ] on both sides of 
Eq.(0 solves the problem with: 



1 kl 1 + ak 2 V s + ik 3 V s 



In the asymptotic limit {kr 3> 1), one obtains: 



%(f ) ~ exp(z/c.f ) + 3 (e k .e r ) fi 



exp(ifcr) 
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+ a + ik. and ik — k/k. (8) 



Eq.JSJl coincides exactly with the general expansion of 
the inverse p-wave scattering amplitude f\ at second or- 
der in the low energy limit (20j , showing that the pseudo- 
potential provides a modeling of two-body collisions 
for relative momenta k <C a. The zero-range potential 
approach requires also that kR <C 1 where R denotes the 
potential range, that is the radius such that for r > R 
the 'real' potential term experienced by particles can be 
neglected in the Schrodinger equation. In actual exper- 
iments on spin-polarized fermions 0, 0, an external 
magnetic field B tunes the energy of a two-body bound 
state in a closed channel which is coupled to the open 
channel associated with the two colliding atoms. A Fesh- 
bach resonance occurs for a vanishing value of the bound 



state energy at a magnetic field B = Bo- Close to 
the resonance, the scattering volume takes arbitrarily 
large values with V s oc —1/(B — Bo), while t he p aram- 
eter a is a slowly varying function of B [|| l25j |. The 
two-body bound state in the closed channel extends on 
short lengths of the order of the potential range and is 
not described by the zero-range approach. However, in 
the regime where a 3 V s is large and positive (B < Bo), 
the pseudo-potential supports a weakly bound state of 
energy es — —h 2 /maV s oc B — Bo which corresponds to 
the outer part (region r > R) of the molecular state 
resulting from the coupling between the closed and the 
open channels. We anticipate that this state is populated 
by pairs of particles in the BEC region of the BEC-BCS 
crossover regime. For large and negative values of a 3 V s 
(B > Bo), this state transforms into a long lived quasi- 
bound state poj . 

Having in hand the short range behavior of the wave 
functions in Eq.JSJ), one may wonder what is the under- 
lying structure of the Hilbert space spanned by the two- 
body eigenstates of the zero-range Hamiltonian. As a 
consequence of the 1 /r 2 singularity, the low energy bound 
state is not normalizable. Moreover, the usual scalar 
product between two different scattering states gives an 
infinite result. Nevertheless, as we show in the following, 
it is possible to introduce a regularized scalar product 
such that the eigenfunctions arc orthogonal to each oth- 
ers and the low energy bound state is normalizable. For 
this purpose, we perform the scalar product between two 
states \^>^) and with k ^ k', but exclude from inte- 

gration the inner volume of the sphere S r , defined above. 
One finds, as r — > 0: 



{r-ar 2 )[ d 2 f'm(f')^ p (r') + 0{r). (9) 
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This integral is non zero and diverges in the zero-range 
limit due to the p-wave singularity of the scattering 
states. Therefore the zero-range Hamiltonian is not her- 
mitian with respect to the usual scalar product. This 
feature follows from the fact that the mapping between 
the true scattering states associated with the finite range 
potential experienced by particles and the states {l^}} 
in Eq.© is not valid for r < R while it is justified out- 
side the potential range. Indeed, the singular boundary 
behavior imposed on wave functions in Eq. @ is the way 
to reproduce the effect of the true finite range potential 
for r > R but has a formal character for interparticle 
distances r < R. In order to conciliate the zero-range 
pseudo-potential with a finite range potential of small 
radius R, we consider the sphere Sr which separates the 
outer part (region r > R) and the inner part (region 
r < R) of the true wave functions. Contribution in the 
scalar product of the outer parts of the true wave func- 
tions is given by Eq.@ with r — R and due to orthogo- 
nality between the true scattering states, the scalar prod- 
uct of the inner parts cancels this term. By construction, 
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the inner parts of the wave functions (r < R) are not de- 
scribed in the zero-range approach. However, we are free 
to modify the usual scalar product in order to take into 
account their contribution. Therefore, we define a regu- 
larized scalar product ( . | . )o by subtracting the surface 
term (r.h.s. of Eq.©) from the usual scalar product. In 
the formal zero-range limit, the regularized scalar prod- 
uct can be also written with a weight g(r) given by: 



g(r) = 1 + 8(r) [(ar 2 — r) . 



(10) 



so that finally, for wave functions obeying the boundary 
condition Eq.®, the scalar product reads: 



(11) 



One can check that with this new metrics, the zero-range 
Hamiltonian is hermitian in the domain defined by wave- 
functions satisfying Eq. @ . Note that following the same 
reasoning, the notion of regularized scalar product has 
been generalized recently for zero-r ang e interactions in 
all the other partial wave channels |26( in the resonant 
regime. Coming back to the low energy bound state of 
energy eg = —h 2 n 2 B /m with n B 2 ~ aV s appearing in the 
regime V s a 3 1, the radial part normalized with respect 
to Eq.lQIj! is: 



n B (r) 



i/a - 3kb/2 



exp(-K B r) 



(12) 



As a result, the probability of finding the molecule of van- 
ishing binding energy outside the range of the true poten- 
tial is of the order of 1/ aR, meaning that the approach is 
consistent for a > and aR > 1 (for example aR ~ 2.8 
in 40 K for the two resonances with B ~ 198.5 G 0]). 
The normalization factor in Eq. (|12fl can be also deduced 
from the residue of the scattering amplitude JSJ) at the 
energy h 2 k 2 /m = eb [2(| and the regularized scalar prod- 
uct extends this method to inhomogeneous situations. 

Now we turn to the case of a homogeneous system 
made of N spin-polarized identical fermions of mass to. 
For a density (n) sufficiently small that binary processes 
occur at low energy (nR 3 < 1) the zero-range approach 
is justified and the interaction between two given parti- 
cles can be modeled by the pseudo-potential in Eq.QJ 
|27j . For a qualitative picture of the physics involved 
in the neighborhood of the resonant regime, we use the 
box model introduced in Ref.|28| in the context of s-wave 
scattering. In this effective approach, the interaction of 
one fermion with all the others is modeled by the inter- 
action of a fictitious particle of mass equal to the reduced 
mass m/2, with a fixed scatterer at the center of a spher- 
ical box of radius L. The wave function of the fictitious 
particle which is non-zero in the p-wave channel only is 
an eigenstate of the pseudo-potential Q and vanishes 
on the surface of the box. It represents the pair func- 
tion of two fermions in their center of mass frame. The 
boundary condition at r = L, mimics the effects of cor- 
relations between pairs and the radius L is fixed by con- 
sidering the non interacting case. In this situation, the 
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FIG. 1: Energy per particle for a homogeneous medium 
composed of polarized fermions as a function of — 1/V S for 
a = lOfcF (dashed line) and a — 10 3 Uf (continuous line). In- 
set: equation of state at resonance (B = Bo). In this regime, 
the energy per particle varies linearly with the atomic density. 



total energy of the gas is E = ^Ncf, where the Fermi 
energy tp = h 2 kp/2m is related to the mean density n, 
through kj? = 6ir 2 n. The gas energy E is also related to 
the energy e of the fictitious particle by E — Ne/2. Fi- 
nally, we obtain a relation between the radius L and the 
Fermi momentum kp by considering the ground state en- 
ergy of the fictitious particle in the box where fcp-L ~ 5.8. 
Hence, L is as expected of the order of the mean inter- 
particle distance. In the interacting case, the equation of 
state is deduced from the energy of the fictitious particle 
e = h 2 k 2 /m = 2E/N, where k (real for e > or imagi- 
nary otherwise) is a solution of the equation: 



kLcos(kL) — sin(fcL) V s k 3 
kL sin(fcL) + cos(fcL) ~ 1 + aV s k 2 



(13) 



In current experiments 000 if we neglect the trap ge- 
ometry, the density and the external magnetic field are 
the two control parameters. We consider then solutions 
of Ea. l|13|) for given values of — 1/V S oc B — Bo and of the 
dimensionless parameter a/kp ^> 1. In Fig.QJ, we have 
plotted the two first branches of the energy per parti- 
cle as a function of — 1/V S for two values of the density 
(a/kF = 10 and 10 3 ). The left part of the upper branch 
corresponds to the metastable weakly repulsive atomic 
phase. The right part of the ground branch is associated 
with the weakly attractive atomic phase where a BCS 
phase is expected at sufficiently low temperature. The 
left part of the ground branch represents the molecular 
phase composed of dimers of energy e/ep ~ —2/V s ak F 
|29| . As explained in Ref.|2^|, this two branches picture 
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provides the two possible scenarios for obtaining a molec- 
ular phase from the weakly interacting atomic Fermi gas 
by varying the external magnetic field. In the first sce- 
nario the system is initially a weakly attractive Fermi 
gas and follows the ground branch while the parameter 
— 1/Vs is tuned from positive to negative values. In the 
second scenario, the system is initially a weakly repul- 
sive Fermi gas and is driven in the resonant regime by 
increasing values of — 1/V S . The ground branch is then 
populated by three body recombination processes. At 
resonance {B = B ) the ground state energy varies lin- 
early with density: E/N oc h 2 n/ma <C cf (see inset of 
Fig-Q). This result which has been found also in a BCS 
treatment strongly differs from the n 2 / 3 law obtained 
in the unitary regime for two-spin components fcrmions 
|28| . Qualitatively, a pair of interacting fermions which 
would have zero energy and infinite size at resonance in 
absence of other particles, is confined by its neighbors in 
a volume of the order of L 3 oc 1/n and the box model 
provides an estimation of the cost in energy of this config- 
uration which is ~ h 2 /maL 3 <C h 2 /mL 2 . Also the level 
crossing between the two branches differs from the s- 
wave result: it is not only translated to the right of the 
zero energy resonance (B = Bo) but also its amplitude 
tends to zero for decreasing densities. As a consequence 



one may expect that in the region of the level crossing, 
non adiabatic transfers of atoms from the atomic branch 
to the ground branch resulting from inelastic three-body 
collisions is followed by only a small heating of the gas 
due to the small kinetic energy of the outgoing resulting 
states. 

To conclude, we comment on the stability of the molec- 
ular phase which is one of the main issue for achieving 
the BEC-BCS crossover. Lifetime of the shallow diatomic 
molecules is limited by recombinations into deep bound 
states of energies 0(—h 2 /mR 2 ). Following, the analy- 
sis of Ref.flTj the corresponding loss rate can be esti- 
mated from the probability that three atoms are con- 
fined within a volume of characteristic length R. The 
fact that the probability for a pair to be in the inner 
part of the potential (r < R) is given by ~ (1 — 1/a.R) is 
a strong indication that this rate depends on the width 
of the resonance considered and could be reduced in the 
cases where aR ~ 1. In contrast, these losses are uni- 
versally suppressed in the unitary two-spin components 
Fermi gas [l7| . 
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